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1 . Introduction 

Besides its simple molecular structure, the magic of 2D graphene, a sheet of carbon graphite, is 
essentially due to two fundamental electronic properties: First for its peculiar band structure 
where valence and conducting bands intersect at two points K+ and K- of the reciprocal 
space of the 2D honeycomb making of graphene a zero gap semi-conductor. Second, for 
the ultra relativistic behavior of the charge carriers near the Fermi level where the energy 
dispersion relation E = E (p) behaves as a linear function in momenta; E (p) = v^p + O (p^) . 
This typical property, which is valid for particles with velocity comparable to the speed of 
light, was completely unexpected in material science and was never su spected before 2004 : 
the year where a sheet of 2D graphene has been experimentally isolated iCeim & Novoselovl. 
I2OO7I; iNovoselov et alJ. l2004h . From this viewpoint, graphene is then a new material with 
exotic properties that could play a basic role in the engineering of electronic devices with 
high performances; it also offers a unique opportunity to explore the interface between 
condensed interface between condensed matter physics and relativistic Dirac theory where 
basic properties like chirality can be tested; and where some specific features, such as 
numerical simulati on methods, can be mapped to 4D lattice gauge theory like lattice QCD 
iBorici, 2008; Capi tani etakl, I2OIOI; ICreutj, l2008h . Although looking an unrealistic matter 
system, interest into the physical properties of graphene has been manifested several decades 
ago. The first model to analyze the band structure of graphite in absence of externa l fields 
was developed by Wallace in 1947 iWallacel Il947l) ; see also iSlonczewski & Weissl |l958). 
Since then, several theoretica l stud i es have been p erformed on graphene in the presence 
of a magnetic field jHaldanel ll983h - jGoerbig et all I2OO6) . The link between the electronic 
properties of graphene and (2 + 1) -dimensional Dirac theory was also considered in many 
occasions; in particular by Semenoff, Fradkin and Haldane durin g the 80-th of the last century 
(Castro-Neto et al., 2009; Haldane, 1988; Semenoff, 1984); see also Ilackiw & Pj 2007: 2008) and 
refs therein. 

In this book chapter, we use the tight binding model as well as the SU (3) hidden symmetry 
of 2D honeycomb to study some physical aspects of 2D graphene with a special focus on 
the electronic properties. We also develop new tools to study some of graphene's cousin 
systems such as the ID- poly-acetylene chain, cumulene, poly-yne, Kekule cycles, the 3D 
diamond and the 4D hyperdiamond models. As another application of the physics in higher 
dimension, we also develop the relation between the so called four dimensional graphene first 
studied in ^Bedaque et al... .2008.: .Creutz.. ,2008.: .Drissi et aL ,2011 a|) : and 4D lattice quantiom 
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chromodynamics (QCD) model considered recentl y in the lattice quantum field theory (QFT) 
literature to deal with QCD numerical simulations iCapitani et al.l.l2009 allb|l . 
The presentation is as follows: In section 2, we review the main lines of the electronic properties 
of 2D graphene and show, amongst others, that they are mainly captured by the SLf(3) 
symmetry of the 2D honeycomb. In section 3, we study higher dimensional graphene type 
systems by using the power of the hidden symmetries of the underlying lattices. In section 

4, we give four examples of graphene's derivatives namely the ID- poly-acetylene chain, 
having a SU{2) invariance, as well as Kekule cycles thought of as a particular ID- system. 
We also study the 3D diamond model which exhibits a Sii (4) symmetry; the corresponding 
2D model, with S!i(3) invariance, is precisely the graphene considered in section 2. In section 

5, we develop the four dimensional graphene model living on the 4D hyperdiamond lattice 
with a SU(5) symmetry. In section 6, we study an application of this method in the framework 
of 4D lattice QCD. Last section is devoted to conclusion and comments. 



2. Two dimensional graphene 

First, we give a brief review on the tight binding modeling the physics of 2D graphene; then 
we study its electronic properties by using hidden symmetri es. We show amongst others that 
the 2D honeycomb is precisely the weight lattice oi SU (3) jPrissi et all l201oh : and the two 
Dirac points are given by the roots of SU (3). This study may be also viewed as a first step 
towards building graphene type systems in diverse dimensions. 

2.1 Tight binding model 

Graphene is a two dimensional matter system of carbon atoms in the sp^ hybridization 
forming a 2D honeycomb lattice. This is a planar system made of two triangular sublattices 
A2 and B2', and constitutes the building block of the layered 3D carbon graphite. Since its 
experimental evidence in 2004, the study of the electronic properties of graphene with and 
without external fields has been a bi g subject of interest; some of its main physical aspects 
were reviewed in iCas tro-Neto et all 12009) and refs therein. This big attention paid to the 
2D graphene, its derivatives and its homologues is because they offer a real alternative for 
silicon b ased technology and bring together issues fr om condensed matter and high energy 
physics iGiuliani et al. , 2010l)jChakrabarti et al.l j2009h allowing a better understanding of the 



electronic band structure as well as their special properties. 
In this section, we focus on a less explored issue of 2D graphene by studying the link 
between specific electronic properties and a class of hidden symmetries of the 2D honeycomb. 
These symmetries allow to get more insight into the transport property of the electronic 
wave modes and may be u sed to approach the defects and the boundaries introduced 
in the graphene monolayer jCortijo & Vozmedianol. l2009h . The existence of these hidden 
symmetries; in particular the remarkable hidden SU (3) invariance considered in this study, 
may be motivated from several views. For instance from the structure of the first nearest 
carbon neighbors like for the typical (Aq-Bi), (A0-B2) / (^0-^3) as depicted in triangle of 
fig{l). These doublets Ag-Bi, AQ-B2, Aq-Bj are basic patterns generating the three S(J(2) 
symmetries contained in the hidden SU (3) invariance of honeycomb. The A-B patterns 
transform in the isospin | representations of SU (2) and describe the electronic wave doublets 
4>± 1 = ['^ (r) / (r)] interpreted as quasi-relativistic 2D spinors in the nearby of the Dirac points 
ICastro-Neto et al],l2009l) . The SU (3) hidden symmetry of honeycomb is also encoded in the 
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Fig. 1. Sublattices A and B of the honeycomb with unit cell given by dashed area. A- type 
carbons are given by red balls and B-type atoms by blue ones. Each carbon has three first 
nearest neighbors as shown by the triangle; and six second nearest ones. 

second nearest neighbors ((Aq-A,)) and ((Bq-B,)) ,i = l, 6 which capture data on its adjoint 
representation where the six ((Ag-A,)) (and similarly for ((Bq-B,))) are precisely associated 
with the six roots of SLf (3) namely iaj, ±0:2, ±0:3; see below. In addition to above mentioned 
properties, hidden symmetries of graphene are also present in the framework of the tight 
binding model with hamiltonian, 

3 / N 

H=-tEL arM+v„ -i' E flr,4 + fcr.fcj, ]+hc , (2.1) 

r,- )j=l r,-,r^ ^ ' ' ' 

where i — 2.8eV is the hopping energy; and where the fermionic creation and annihilation 
operators a, b, a' , b'^ are respectively associated to the pi-electrons of each atom of the 

sublattices A2 and B2. The three relative vectors vj, V2, V3 define the first nearest neighbors, 

see fig© for illustration. These 2D vectors are globally defined on the honeycomb and obey 

the remarkable constraint equation 

vj + V2 + V3 = 0, (2.2) 

which, a priori, encodes also information on the electronic properties of graphene. 
Throughout this study, we show amongst others, that the three above mentioned SU (2)'s 

are intimately related with these v,, 's which, as we will see, are nothing but the weight vectors 

A„ of the SU (3) symmetry; i.e v„ = ''^x^- The wave functions cp^^^^ (r) of the delocalized 

electrons are organized into a complex SU (3) triplet of waves as given below 

/|Ai)\ 

IA2) ^3, A1 + A2 + A3 =0. (2.3) 
\|A3>/ 

The symbol 3 refers to the 3-dimensional representation of S !i (3); say with dominant weight 
Aj. We also show that the mapping of the condition E;,^]^ A„ = to the momentum space can 
be interpreted as a condition on the conservation of total momenta at each site of honeycomb. 
This connection with SU (3) representations opens a window for more insight into the study 
of the electronic correlations in 2D graphene and its cousin systems by using symmetries. 
The organization of this section is as follows: In subsection 1, we exhibit the SU (3) symmetry 
of graphene. We also give a field theoretic interpretation of the geometric constraint equation 
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vi + V2 + V3 = both in real and reciprocal honeycomb. We also use the simple roots and the 
fundamental weights of hidden SU (3) symmetry to study aspects of the electronic properties 
of 2D graphene. In subsection 3, we develop the relation between the energy dispersion relation 
E [kx,ky) and the hidden SU (3) symmetry. Comments regarding the link between graphene 
bilayers and symmetries are also given. 

2.2 Symmetries and electronic properties 
2.2.1 IHidden symmetries of graphene 

In dealing with pristine 2D graphene, one immediately notices the existence of a hidden 
S(J (3) group symmetry underlying the crystallographic structure of the honeycomb lattice 
and governing the hopping of the pi-electrons between the closed neighboring carbons. To 
exhibit this hidden SU (3) symmetry, let us start by examining some remarkable features 
on the graphene lattice and show how they are closely related to S(J(3). Refereing to the 
two sublattices of the graphene monolayer by the usual letters A2 and B2 generated by the 
vectors ai = d(\/3,0), a2 = j{-\/3,3); together with the three relative vi = |(\/3,1), 
V2 = V^S, 1), V3 = — vj — V2 with carbon-carbon distance d ~ 1.42 A; and denoting by 



(py^ (r,) and cpg yijj the wave functions of the corresponding pi-electrons, one notes that the 
interactions between the first nearest atoms involve two kinds of trivalent vertices capturing 
data on S(J(3) symmetry, see figlE) for illustration. This hidden SU (3) invariance can be 



Fig. 2. (a) Nearest neighbors of a A- type atom, (b) Nearest neighbors of a B-type atom. 
These two configurations are precisely the representations 3 and 3* of SU (3) . 

made more explicit by remarking that the relative vectors vj, V2, V3 describing the three first 
closed neighbors to a A- type carbon at site r, of the honeycomb, together with their opposites 
— v„ for B-type carbons, are precisely the weight vectors of the 3-dimensional representations 



of the S(J(3) symmetry, v„ = dJj\„. For readers not familiar with representation group 



theory terminology, we give here below as well in the beginning of section 3. a summary on 
the SU (N) symmetry. 
some useful tools on SU{3) 

Roughly, the SU (3) symmetry is the simplest extension of the SU (2) symmetry group behind 
the spin of the electron. The basic relation vi + V2 + V3 = of the honeycomb, which 

upon setting v„ = d^J^\„, reads also as Ai + A2 + A3 = 0. This constraint relation has an 
interpretation in SLf (3) representation theory; it should be put in one to one correspondence 
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with the well known SU (2) relation — j) = of the spin j representation; 

-1^ 



1-1 = 0, (2.4) 



see also eql l3.3t for details. The basic properties of the SU (3) symmetry are encoded in the so 
called Cartan matrix and its inverse KJ^^ which read as 



These matrices can be also written as the intersection of 2D- vectors as Kjj = a; ■ aj, K^^^ = 
(Vi ■ Wj where and cc2 are the two simple roots of Sfi (3) and where wi and iV2 are the 
corresponding two fundamental weights which related to the simple roots by the following 
duality relation 

Oii ■ Wj = Sij , a, = -KyWy. (2.6) 

Using these tools, the honeycomb relation Ai + A2 + A3 = is naturally solved in terms of the 
fundamental weights as follows 

Ai = A2 = ix'2 — (jJ\, A3 = 61^2 . (2.7) 

We also have the following relations between the a, vectors and the v; ones: aj = (vj — V2), 
32 = V2 — V3 and 33 = V3 — vj. Notice that the vectors ±ai, ±a2, ±33 are, up to the scale 

factor A\J\, precisely the six roots of the SLf (3) symmetry 

'3 



ai = d\^\cc\, 32 = 2^2, 33 = — f (ai + 0C2) , (2-8) 

where we have also used the remarkable relation between roots and weights that follow from 
eqs <2.5l2.6t . 

ci\ = liVi — LO2, CX.2 = 2oj2 — (2.9) 



2.2.2 Electronic properties 

Quantum mechanically, there are two approaches to deal with the geometrical constraint 
relation ( 12.21 . The first one is to work in real space and think about it as the conservation 
law of total space-time probability current densities at each site r, of the honeycomb. The 
second approach relies on moving to the reciprocal space where this constraint relation 
and the induced electronic properties get a remarkable interpretation in terms of Sii(3) 
representations. 

1) conservation of total current density 

In the real space, the way we interpret eq (|2.21 is in terms of the relation between the time 
variation of the probability density p (t, r,) = (t, r,) | of the electron at site r, and the sum 
I^H=i Jv„ (','■/) = J (f, >■(■) of incoming and outgoing probability current densities along the v„- 
directions. On one hand, because of the equiprobability in hopping from the carbon at r,- to each 
one of the three nearest carbons at r, + v,,, the norm of the Jv„- vector current densities should 
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be equal and so they should have the form 

Jv„(f,r,)=/(t,r,)e„,« = l,2,3. (2.10) 

These probability current densities together with the unit vectors e„ = ^ pointing in the 
different v„- direction; but have the same non zero norm: ||Jvi|| = ||Jv2l| = IIJV3II = |/|- 
Substituting in the above relation, the total probability current density J {t, r,) at the site r and 
time t takes then the factorized form 




(2.11) 



On the other hand, by using the Schrodinger equation ih^ = ^ + V j 1^ describing the 
interacting dynamics of the electronic wave at r, we have the usual conservation equation, 

^^|^+divj(f,r) =0 , (2.12) 

with probability density p {t, r) as before and ] = {^Vtp* — (p*V(p) with m the mass of the 
electron and (p = cp {t,i) its wave. Moreover, assuming ^ = corresponding to stationary 
electronic waves (p {t,i) = e'"'(p (r), it follows that the space divergence of the total current 
density vanishes identically; divj = 0. This constraint equation shows that generally J 
should be a curl vector; but physical consideration indicates that we must have J (^, r) =0, in 
agreement with Gauss-Stokes theorem Jy div J dV = J^y J.dcr leading to the same conclusion. 

Combining the property J (t, r) =0 with its factorized expression ^ (X]„ v„) given by eq O.llt 
together with j 7^ 0, we end with the constraint relation v„ = 0. 
2) conservation of total phase 

In the dual space of the electronic wave of graphene, the constraint relation <2.2t may be 
interpreted in two different, but equivalent, ways; first in terms of the conservation of the 
total relative phase A<pj„j = X^k.Ar of the electronic waves induced by the hopping to the 
nearest neighbors. The second way is in terms of the conservation of the total momenta at 
each site of the honeycomb. 

Decomposing the wave function tp (r), associated with a A-type carbon at site r, in Fourier 
modes as J^i^e'^"^'' ^ (k); and similarly for the B-type neighboring ones i^(r + v„) = 
e'-^"^^ '^ ^„ (k) with k = (fcx, /cy) , we see that f (k) and the three (fc) are related as 

(fc) = e'^'^""^ (k) , H = 1,2,3 , (2.13) 

with relative phases 6,-, = k • v„. These electronic waves have the same module, (k) |^ = 
1^ (fc)l^; but in general non zero phases; 9\ ^ di ^3- This means that in the hop of an 
electron with momentum p = ?!k from a site r, to the nearest one at r; + v„, the electronic 
wave acquires an extra phase of an amount 0„; but the probability density at each site is 
invariant. Demanding the total relative phase to obey the natural condition. 



01 + 02 + 03 = , mod (27r) , 



(2.14) 
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one ends with the constraint eq l|2.2> . Let us study two remarkable consequences of this special 
conservation law on the 9n phases by help of the hidden S (J (3) symmetry of graphene. Using 
eq d2.7t , which identifies the relatives v„ vectors with the weight vectors Wn, as well as the 
duality relation «, ■ Wj = Sjj (2|6l, we can invert the three equations 6,, = k ■ v„ to get the 
momenta p„=?ik,; of the electronic waves along the v„-directions. For the two first 9„'s, that 
is n = 1, 2, the inverted relations are nicely obtained by decomposing the 2D wave vector k 
along the cci and 1^2 directions; that is k = fciai + ^2*^2; arid end with the following particular 
solution, 

01 = kid, 02 = {k2 - fci) d, 63 = -k2d . (2.15) 



2.3 Band structure 

We first study the case of graphene monolayer; then we extend the result to the case of 
graphene bilayers by using the corresponding hidden symmetries. 

2.3.1 Graphene monolayer 

By considering a graphene sheet and restricting the tight binding hamiltonian J2.lt to the first 
nearest neighbor interactions namely, 

2 

H = -t'£J^ar,bl+,_,+hc, (2.16) 

r, n=l 

we can determine the energy dispersion relation and the delocalized electrons by using 

the S[i(3) symmetry of the 2D honeycomb. Indeed performing the Fourier transform of 
the various wave functions, we end with the following expression of the hamiltonian in the 
reciprocal space 




(2.17) 



The diagonalization of this hamiltonian leads to the two eigenvalue E± = ±t |£k| giving the 

energy of the valence and conducting bands. In these relations, the complex number is an 
oscillating wave vector dependent function given by = e"^^^ + e"^'^'^ + e^"^^Q'^+Q^-) where 
we have set Q/ = k.A;. This relation, which is symmetric under permutation of the three Q,, 
can be also rewritten by using the fundamental weights as follows, 

= ^e'rf[k.Wi] ^ f,id\k.(w2-wi)] _|_g-/rf[k.a;2]j ^ (2.18) 

Up on expanding the wave vector as k = ki^i + k2K2, this relation reads also as = 
gikid _|_ ^i(ki-h)d _|_ g-fizd Notice that from l l2.18t , we learn that is invariant under the 
translations k — > k+ ^ (Njaj + N2CC2) with Nj, N2 arbitrary integers; thanks to the duality 
relation Xi.LOj = Sjj. 

Notice also that near the origin k = 0, we have = 3 + O (k^) , in agreement with non 
relativistic quantum mechanics. The three terms which are linear terms in k cancel each 
others due to the S(J (3) symmetry. Notice moreover that the Hamiltonian l |2.16t has Dirac 
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Fig. 3. On left the band structure of the graphene monolayer where one recognizes the Dirac 
points. On right, it is shown the relativistic behavior near a Dirac point where conducting 
and valence bands touch. On right the band structure in GGA approximation using QE 
code. 

zeros located, up to lattice translations, at the following wave vectors 



Kf = ±|f {Iwi - W2) , = ±|f (-wi + , = ±|f (wi + a;2) , (2.20) 



are not completely independent; some of them are related under lattice translations. For 
instance, the three are related to each others as follows 



These two zeros are not related by lattice translations; but are related by a Z2 symmetry 
mapping the fundamental weights and the simple roots to their opposites. 
We end this section by noting that the group theoretical approach developed in this study 
may be also used to deal with graphene multi-layers and cousin systems. Below, we describe 
briefly the bilayers; the cousin systems are studied in next sections. 

2.3.2 Bilayer graphene 

Bilayer graphene was studied for the first time in lMcCann & Falkol j2006|) . It was modeled 
as two coupled hexagonal lattices including inequivalent sites in the two different layers that 
are ranged in the Bernal stacking {ihe stacking fashion of graphite where the upper layer has its B 
siiblattice on top of sublattice A of the underlying layer) as showed in the figure (S). This leads 
to a break of the Dg;, Bravais symmetry of the lattice with respect to the c axis. Comparing 
bilayer graphene to monolayer one, we notice that its unit cell contains four atoms. There 
exist other arrangements such as the AA stacking, where the two lattices are directly above 
each other and bonds form between the same sublattices. The AB stacking arrangement 




(2.19) 



Notice that these six zero modes, which read also as 



K+ + ^Ju,2 = K+ + f 0.1 = K3+ . 



(2.21) 
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Fig. 4. On right bilayer graptiene in the AB stacking allotrope with hop energy t^ between 
the layers. On right bilayer graphene projection in the x-y plane. 



was experimentally verified in epitaxial graphene by Ohta et al. jOhta et all |2007|) . The 
tight-binding model describing bilayer graphene is an extension of the one corresponding to 
the monolayer {2.H , by adding interlayer hopping elements H = H-i+ H2 + where H, are 
as in il.H and where 



H± = -t^ZLai (r,) fo| (r;) + ^2 (r,) (r,) + he , 

n=l 



(2.23) 



.with f I is the hon enerey of the pi-electrons between layers calculated to be f 1 — Tn 

iCharlie^r et alJ.Tl99lh . - -■ - ■' - ' " ' - - ' - - 



From the view of hidden symmetries, the bilayer graphene has a 
symmetry type SU (3) x S!i(2) x SLf (3); each Sfi(3) factor is associated with a graphene 
sheet; while the SU (2) corresponds to the transitions between the two layers and is associated 
with propagation along the z-direction of the 3D-space. 

Applying Fourier transform, the above hamiltonian can be rewritten in the following form: 



k 



/ 











f flik\ 







t± 

t 












t± 
t 









«2k 


V-f 







J 




\fc2k/ 



(2.24) 



with is as in eq( l2.18t . The diagonalization of this hamiltonian leads to the following energy 
dispersion relations. 



(2.25) 



The corresponding band structure has two additional bands, n and n* states having 
lower energy bands, that i s consequence of the n umber of atoms per unit cell. Neutral 

bilayer graphene is gapless iMcCann & F^ b006h and exhibits a variety of second-order 

effects. The studies on bilayer graphene show that it has many common physical properties 

with the monolayer, such as the exceptionally hig h electron mobility and high mechanical 

stability jOhta et all, l2007l) - lN"ovoselov et alj j2006l) . The synthesis of bila yer graphe ne thin 

films was realized by deposition on a silicon carbide (SiC) substrate l lOhta et all 12006). 

The measurements of their electronic band structure, using angle-resolved photo-emission 

spectroscopy (ARPES), suggest the control of the gap at the K Point by applying Coulomb 

potential between the two layers. This tuning of the band gap changed the biased bilayer 

from a conductor to a semiconductor. 
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r KM r 

Fig. 5. Band structure of the bilayer graphene 

3. Higher dimensional graphene systems 

Motivated by the connection between 2D graphene and SU (3) symmetry, we study in this 
section the extension of the physics of 2D graphene in diverse dimensions; that is ID, 2D, 3D, 
4D, and so on; the 2D case is obviously given by 2D graphene and its multi-layers considered 
in previous section. The precited dimensions are not all of them realizable in condensed matter 
physics; but their understanding may help to get more insight on the specific properties of 2D 
graphene since the SU (3) is the second element of the SU (N) symmetries series. 
First we develop our proposal regarding higher dimensional graphene systems that are based 
on SU (N) symmetry including the particular ID poly-acetylene chain which corresponds to 
SU (2) symmetry. Then, we compute the energy dispersion relation of these kinds of lattice 
quantum field theory (QFT). Explicit examples of such lattice fermionic models will be studied 
in the next sections. 

3.1 The Sfi(N) model 

Higher dimensional graphene systems are abstract extensions of 2D graphene; the analogue of 
the 2D honeycomb is given by a real N-dimensional lattice Cc^^^i^ f^y The quantum hamiltonian 
describing these systems is a generalization of (|2.11 and reads as follows, 

Hn = -tE f E ar.bl+A +hc-t'Z E («r,4.+v„,„ + ^r^l+Y ) ' (3-1) 

r,- \«=1 / r; n<m=l ^ ' 

where «r,v i'i,+v„/ Hi+v„ fermionic annihilation and creation operators living on -Cs„(jv)- 
Moreover the vectors vj, v^r are, up to a global scale factor, the fimdamental weights of the 
N-dimensional representation of the SU (N) symmetry constrained by the typical property 

vi + ... + VN = 0. (3.2) 

The vectors Vnm = (vn — v,„) are, up to a scale factor, precisely the N (N — 1) roots of SU (N); 
they obey as well the group property E = 0. 

These particular features of Hjv let understand that its physical properties are expected to be 
completely encoded by the hidden SU (N) symmetry of the model. Below, we show that this 
is indeed the case; but for simplicity we will focus on the first term of Hfj; i.e working in the 
limit t' 0. 
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3.1.1 Useful tools on SU (N) symmetry 

Since one of our objectives in this paper is to use the S(J(N) symmetry of the crystals to 
study higher dimensional graphene systems; and seen that readers might not be familiar with 
these tools; we propose to give in this subsection some basic tools on SLf (N) by using explicit 
examples. 

a) cases SU {2) and S(J (3) 

The S U (2) symmetry is very familiar in quantum mechanics; it is the symmetry that describes 
the spin of the electrons and the quantum angular momentum states. 

Roughly speaking, the SLf (2) symmetry is a 3-dimensional space generated by three matrices 
which can be thought of as the usual traceless Pauli matrices 

involving one diagonal matrix o"", giving the charge operator, and two nilpotent matrices 

interpreted as the step operators or equivalently the creation and annihilation operators in the 

language of quantum mechanics. These three matrices obey commutation relations 

±2a^ that define the sm (2) algebra. Observe also the traceless property of the charge operator 

Trcr" = ^ — ^ = 0, which should be related to the constraint relation l l3.2t with N = 2. 

The SU (3) symmetry group is 8-dimensional space generated by 8 matrices which can be 

denoted as 

^1,^2, , (3.4) 

with hi, /i2 two diagonal matrices defining the charge operators and six step operators e*"i, 
g±a2^ giffa piayij^g jole of creation and annihilation operators. The e*"''s are nilpotent and 
are related as (e^"' )^ = . An example of these matrices is given by the following 3x3 
matrices 



(3.5) 

















[l 


f'2 




, h2=[ 




») 







F3/ 


\ 











with the traceless property of the charge operators which reads as follows 



Ai + A2 + A3 = , A,- = 




and which should be compared with the case N = 3 in <3.2t . 

The vectors «i and «2 are the simple roots encountered in the previous section; their scalar 
product ccj.ccj gives precisely the Cartan matrix K,y of eq <2.5t . 

b)caseSU{N) 

In the general case N > 2, the corresponding SU{N) symmetry is (N^ — 1) -dimensional 
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space generated by (N^ — l) matrices; N — 1 of them are diagonal 

hi, ... 



(3.7) 



and are interpreted as the charge operators; and N (N — 1) step operators giving the creation 
and annihilation operators e^" with a standing for generic roots containing the two 
following: 

(a) N-1 simple ones namely «i,a2,...,ttN-i (together with their opposites) whose scalar 
products ccj.ccj give precisely the following (N — 1) x (N — 1) Cartan matrix 



K 



/ 2 -1 
-1 2 -1 
0-12 




\ 



\ 






2 -1 
-12/ 



(3.8) 



(b) non simple roots given by linear (positive and negative) combinations of the simple ones; 
these roots are given by ±/?,y = ± (^a, ■ + 
Notice that the above Cartan matrix K and its inverse 

N N-1 N-2 



+ Cij ) with 1 < ! < j < N - 1. 



K 



I N+1 N+1 N+1 

N-1 2(N-1) 2(N-2) 

N+1 N+1 N+1 

N-2 2(N-2) 3(N-2) 

N+1 N+1 N+1 



2 4 6 

N+1 N+1 N+1 
\ 1 2 3 

\ N+1 N+1 N+1 



N+1 N+1 \ 



4 

N+1 
6 

N+1 



2 

N+1 
3 

N+1 



2(N-1) N-1 

N+1 N+1 

N-1 N / 

N+1 N+1 / 



(3.' 



capture many data on the SU (N) symmetry; they give in particular the expression of the 
simple roots ai, X2, a^-i in terms of the fundamental weights wi, w^-i and vice versa; 
that is a, = X]/ ^ij^] and w, = ^ij^"'-/- Recall that simple roots and fundamental weights 

obey the duality property «;.a7y = ^;y; we also have Wj.Wj = K^^^. 



3.1.2 The lattice jCsu{n) 

The lattice -C^u^n) a real (N — 1)- dimensional crystal with two superposed integral 
sublattices An and B^; each site fm of these sublattices is generated by the SU (N) simple 
roots ct\, ajv-i; 

rm = rtiioci + m20i2 + ...mjv_iajv-i , (3.10) 

with m,- integers; for illustration see the schema (a), (b), (c) of the figure l|6) corresponding 
respectively to N = 2,3,4; and which may be put in one to one with the sp^, sp^ and sp^ 
hybridization of the carbon atom orbital 2s and 2p. 

On each lattice site r„; of £s;,(iv); say of A-type, lives a quantum state >lr,„ coupled to the 
nearest neighbor states; in particular the first nearest states Br,„+v, and the second nearest 
ones A, +v,- 
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Fig. 6. (a) 1A+2B lattice sites of Csu{2)' A-type in blue and B-type in red; the 2B form a st; (2) 
doublet, (b) 1A+3B sites of J~'sii(3)' ^he 3B form a sii (3) triplet, (c) 1A+4B sites of J~'su{A) with 
4B sites forming a regular tetrahedron. 

Generally, generic sites in Cgn^j^-j have the following properties: 

(1) N first nearest neighbors with relative position vectors v,- constrained as 

vi + . . . + VN = . (3.11) 

These constraint relations are solved in terms of the SU (N) weight vectors A, (resp. —A,) of 
the fundamental (anti-fundamental) representation as follows 

V, = aX, = , (3.12) 

where d is the relative distance between the closest llsu{i>i) sites. The A,'s which satisfy Ai + 
. . . + An = can be nicely expressed in terms of the fundamental weights a;i,...,cfjv_i as 
follows 

Ai = w\, A; = Wi — An = — a;iv-i . (3.13) 

From the QFT view, this means that the quantum states at r„, + v, sites are labeled by the A, 
weights as Br„,+v, = (Th/) and so the multiplet 

/ |Ai > \ 

: ^ N, Ai + . . . + An = 0, (3.14) 

\|An>/ 

transform in the fundamental representation of SU (N). 

(2) N (N — 1) second nearest neighbors of A-type with relative position vectors V,y given by 
v; — vjand obeying the constraint relation Yji,] ^ i] = 0. This condition is naturally solved 
by <3.11t and <3.12t showing that the relative vectors between second nearest neighbors are 
proportional to S!i (N) roots /5,y like 

V,; = , ^ = A, - \j, (3.15) 

and so the condition X] V,y = turns to a SLf (N) property on its adjoint representation labeled 
by the roots. 
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3.2 Energy dispersion relation 

Restricting the analysis to the first nearest neighbors described by eql l3.H in the limit t' — ^ 0, 
the hamiltonian Hj^ on /)sw(N) reduces to 

= -f E f E flr,fcj+v„) + he , (3.16) 

r, Vn=l / 

where now r, and v„ are (N — 1)- dimensional vectors. By using the Fourier transform of the 
field operators flr, and namely, 

fl,, ~ ■ - Ze'^-^''''+^'K (3 17) 

k k 

we can put the hamiltonian Hfj as a sum over the wave vectors k in the following way; 



(3.18) 



with £k = tj^ie'"^-^'. This complex number can be also written as tj^ie'"^-^'^^^^'-'^') with 
cv_i = = wjv. The energy dispersion relation of the "valence" and "conducting" bands are 
obtained by diagonalizing the hamiltonian Hj^; they are given by ± | with, 



l^kl = t 



N 

N + 2 cos [flk. (^A,-AyjJ . (3.19) 

i<j=i 



Notice that |£k| depends remarkably in the difference of the weights A,— A^; which by help of 
eq <3.13t , can be completely expressed in terms of the fundamental weights. 
To get the Fermi wave vectors kp for which the oscillating multi-variable function = 
f e"^^-^' vanish, we will proceed as follows: First, we work out an explicit example; then 
we give the general result. To that purpose, we expand the wave vector k in the weight vector 
basis as follows, 

k = E Q'^-' (Qi.-vQn) eR^, (3.20) 
1=1 

and focus on working out the solution for the particular case where all the Q,'s are equal, 
i.e: Qi = Qi = ■■■ — Qn-1 = Q- General solutions are obtained from this particular case 
by performing lattice translations along the ^/-directions; this leads to the new values Q; = 
Q + with «; integers. Obviously, one may also expand the wave vector k like 

N-l 

k = E ' iki,-,kN) e R~ . (3.21) 

! = 1 

But this is equivalent to <3.20t ; the relation between the Q;'s and the fc;'s is obtained by 
substituting a, = Ey ^ij'^j info j3.211 and identifying it with l |3.20t . To compute the factors 
gidk.A,^ we express the vectors A; in terms of the simple roots as follows 

Aj = loi, \2 = co\ — (Xi, . . . \n = toi — (Xi — ... — (3.22) 
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then we use the root/weight duality relation a7;.ay = Sjj as well as the simple choice 
Q; = Q to put the scalar product k.A; into the following form (k.A/) = (k.wi) — IQ, 
I = 1, N — 1. Putting this expression back into and setting ^ = e'^Q, we obtain 

£k = e'"^-'^^ + ^ + ... + f '^^"^j = which is exactly solved by the N-th roots of unity namely 

Q = ±W'S = 1 ] . (3.23) 

Therefore the Dirac points are, up to lattice translations, located at the wave vectors ks = 

4. Leading models 

In this section, we study the cases N = 2, 4 as N = 3 corresponds precisely to the 2D graphene 
before. The case N = 5 will be studied in the next section seen its remarkable relation with 4D 
lattice QCD. 

4.1 The su (2) model 

In this case, the lattice /^st,(2)/ which is depicted in the figure |7), is a one dimensional chain 
with coordinate positions = ma where a is the site spacing and m an arbitrary integer. 



Aj Aj Aj 




Bj B3 B3 83 B3 



Fig. 7. the lattice £5 „p) given by the superposition of two sublattices .4s„(2) (in blue) and 
^su{2) (ir^ red). The atoms may be thought of as carbons in the sp^ hybridization state. 

Examples of carbon chains with delocalized electrons are given by one of the three following 
molecules 



chain 


molecule 


delocalized electrons 


polyacetylene 




-CH=CH-CH 


= CH-CH- ... 


1 


cumulene 




=c=c=c=c 


= C = ... 


2 


poly-yne 




-C=C-C=C- 


-C = C-... 


2 



(4.1) 



These molecules can be taken as th e graphene bridge ultimately narrowed down to a few- 
carbon atoms or a single-atom width < Giritet et alll2009l:lTurll2008l;lKoskinen et aUboOSl) . Each 
site of J~'su{2) has two first nearest neighbors forming an su (2) doublet; and two second nearest 
ones that are associated with the two roots ±a of su (2) in agreement with the generic result 
summarized in the table. 



nearest neighbors 


SU{N) 


S!i(2) 


SU{3) 


S!J(4) SU{5) 


first 


N 


2 


3 


4 5 


second 


N(N-l) 


2 


6 


12 20 
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In the SU (2) lattice model, eqs <3.21 read as 



vo + vj = , (a) 
Vol = vo - vi , (b) 



(4.3) 



-J of the S!i (2) fundamental 



and are solved by the fundamental weights Aj 
representation; i.e the isodoublet. 
1) poly acetylene 

The hamiltonian of the polyacetylene, where each carbon has one delocalized electron, is given 
by 

= -f £ (flr,>+,+« + «r:»^i-fl) + • (4.4) 



Substituting N = 2 in <3.19t , we get the following energy dispersion relation 



e^l = fi/2 + 2cos(2flfc) 



(4.5) 



which is also equal to It cos {ka) in agreement with the expression e/^ = t (e^"^ + e '"^^ ; see 
also figure Moreover, the vanishing condition £2 (k) = is solved by the wave vectors 
fc± = ±^ mod (^"'^ 




Fig. 8. Energy dispersion relation of ID poly-acetylene chain. 
2) cumulene and poly-yne 

In the case of cumulene and poly-yne, the two delocalized electrons are described by two wave 
functions (p^ ,<p^ . The tight binding hamiltonian modeling the hopping of these electrons is a 
generalization of H(. Let a, , a = 1, 2 (resp. br,„±a' ^t^±a) annihilation and creation 

operators at the site r,,, (resp. r,„ ± a), the hamiltonian reads as follows 



Ht u — 



2 

E E 

m a,/J=l 



(4.6) 
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where f^./^ ^'ci^ hop energy matrices which are identical for cumulene (f^/? — ^^^)/ but 
different for poly-yne {1^^ ^ t'^a). Mapping this hamiltonian to the reciprocal space, we get 



/ 


A 





B\ 




\ 


A* 





c* 










C 





D 


\ 




\B* 





D* 


0/ 




I 



(4.7) 



with 



Aik) = 


hie 




-ika 


B{k) = 




+ V" 


- ika 


C{k) = 






-ika 


D{k) = 


t22e''"' 


+ t'22<^- 


-ika 



(4.i 



Now, using the fact that the two delocalized electrons are indistinguishable, it is natural to 
assume the following relations on the hop energies fn = t22, hi = ill and the same thing 
for the t'^p matrix. This leads to the relations A = D, B = C and so the above hamiltonian 
simplifies. In this case, the four energy eigenvalues are given by 



E± = ± ^{A*+B*) {A + B) 
E'^ = ±y'{A* - B*) {A-B) 



(4.9) 



lika 



and the zeros modes are given by e 

t^a = t' „, it follows that the zero modes are located as A: = ±^ mod ^ 



Since in the case of cumulene we have 



3) nanoruban 

We end this paragraph noting that such analysis may be also extended to the particular case of 
the periodic chain made by the jimction of hexagonal cycles as depicted in the figure This 




Fig. 9. a periodic chain in 3D space with unit cells given by hexagonal cycles. Each cycle has 
six delocalized electrons. 

chain, which can be also interpreted as the smallest graphene nanoruban, is very particular 
from several issues; first its unit cells can be taken as given by the hexagonal cycles; second 
amongst the 6 carbons of the unit cycle, 4 of them have two first nearest neighbors and the 2 
others have three first nearest ones. The third particularity is that the tight binding description 
of this chain is somehow more complicated with respect to the previous examples. Below we 
focus on the electronic properties of a given cycle by using the same approach we have been 
considering in this study. 
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4.2 Kekule cycles 

Kekule cycles are organic molecules named in honor to the German chemist Friedrich Kekule 
known for his works on tetravalent structure of carbon and the cyclic structure of benzene 
QHg. These molecules; in particular the family C„H„ with n > 3, may be thought of as one 
dimensional cycles living in the 3D space; they involve carbon atoms (eventually other atoms 
such as Nitrogen) arranged in a cyclic lattice with both a- and /r-bonds. All these carbon atoms 
are in the sp^ hybridization; they have 3n covalent u-bonds defining a quasi-planar skeleton; 
and n delocalized /r-bonds with Pi electron orbital expanding in the normal direction as shown 
in the examples of figlllOt. Our interest into Kekule molecules, in particular to the C2NH2N 




Fig. 10. Six examples of Kekule cycles type C„ H„ with w = 3, 4, 5, 6, 7, 8 . The cations C+ of 
these molecules form a heavy skeleton represented by n-polygons. The orbitals in the normal 
direction are associated with the delocalized Pi-electrons. 

family, comes from the fact that they can be viewed as the ID analogue of the 2D graphene 
monolayer; they may be also obtained from the poly-acetylene chain by gluing the ends. It is 
then interesting to explore the electronic properties of this special class of systems by using 
the tight binding model and symmetries. To illustrate the method, we focus on the benzene 
CgHg thought of as the superposition of two C3H3 sub-molecules as depicted in figure till . 
From group theory view, the positions vi, V2, V3, V4, V5, V5 of the carbon atoms are given by 





B 



Fig. 11. Kekule molecules as the superposition of two sublattices. Sublattice A in blue and 
sublattice B in red. Except the benzene, these molecules are generally are non planar. 

the six roots of the SU (3) symmetry; that is 



V,- = aoci , V3+,- = -««;, 



(4.10) 
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where a = 1.39 A and where the three a/'s are as in section 2. 
tight binding description 

The electronic properties of the CgHg are captured by the pi-electrons of the carbons. Denoting 
by aj., atj (resp. bj., brj the usual electronic creation and annihilation operators associated 
with the A, (By) atoms in the sublattice Ahem {^benz)> the tight binding hamiltonian of the 
benzene, restricted to first nearest neighbors, reads as follows. 



Hhenz =-t L L \ L ar„„b+ „+v„ ] + hc (4.11) 



3/2 
m=-oo;=i \y=i 

In this relation, the position vectors r;„, have two indices; / and m. The first one takes the 
values I = 1,2,3; it indexes the three atoms in Ahemi ^rid the three ones in Bi,„,z- These 
positions are as follows, 

-A _ -A -A -B _ _B /a i on 

'^2l-\,m ~ '^Im' '^3m' "^Sm ' '^2l,m ~ '^2m' Mm' '^bm ■ \*-^^) 

The second integer is an arbitrary number (m 6 Z); it captures the periodicity of the cycle and 
encodes in some sense the rotational invariance with respect to the axis of the planar molecule. 
To fix the ideas, think about r;,„ as the I-th electron in the sublattice Ai,enz' fhat is r;,„ = r^(_]^ 
After a hop of this electron to the two first nearest carbons in B(,(,„2, the new position is 

4,m=^Un+^lj J = ± ■ (4.13) 

where the v;ys are the relative positions of the first nearest neighbors. 

Taking the Fourier transform of the creation and annihilation operators, = Ejc£*''' '"cJ 
with standing for a^,b^, we get an expression involving the product of three sums X]k Ek' 
Then, using the discrete rotational invariance with respect to the axis of the molecule, 
we can eliminate the sum X],,, in terms of a Dirac delta function S2 (k — k') and end, after 
integration with respect k', with the following result, 

nbcm = L («k, fek) ( °* ( bi ) ■ ^^-^^^ 

with 

3 / N 

£k = - t E (e'^ "'- +e'^ '"^) . (4.15) 

Notice that like in graphene, the above hamiltonian has two eigenvalues ± |£kl. Moreover, 
substituting the v;±'s by their explicit expressions in terms of the SLf (3) roots a/, we obtain 
the following dispersion relation together with a constraint relation capturing the planarity 
property of the molecule 

Ek = -2f E cos f ^^^k.a,) , ;/k = E sin f^^k.a,) = . (4.16) 

;=o ^ ' 1=0 ^ ' 

Notice that the constraint equation y]\^ = Q allows us to express the component of the wave 
vector in terms of k\ and vice versa as depicted in fig<12t. This relation plays a crucial role in 
the determination of the wave vectors at the Fermi level. 
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Fig. 12. the plot of the energy dispersion relation Ei-j ^.^ and the constraint relation j.^ = 
sinfci +sinfc2 — sin(fci+fc2) =0 in the reciprocal space. 

4.3 The diamond model 

The diamond model lives on the lattice J~-su(4)', this is a 3-dimensional crystal given by the 
superposition of two isomorphic sublattices A4 and S4 along the same logic as in the case of 
the 2D honeycomb. Each site r„, in /^s„(4) has 4 first nearest neighbors at (r,„ + v, ) forming the 
vertices of a regular tetrahedron. A way to parameterize the relative positions v, with respect 
to the central position at r,„ is to embed the tetrahedron inside a cube; in this case we have: 



1,-1, +1) , V2 



_d_ 



V3 = ^ (+1,-1,-1) , vo = (+i,+i,+i 



(4.17) 



V3' 



Clearly these vectors satisfy the constraint relation vg + + V2 + V3 = 0. Having these 
expressions, we can also build the explicit positions of the 12 second nearest neighbors; these 
are given by V,y = v, — vy; but are completely generated by the following basis vectors 



Ri = ^ (2,2,0) , R2 = ^ (0,-2,2) , R3 = ^ (-2,2,0) 



that are related to V,y as R, = "V^j^iy 
• the intersection matrix of the R, vectors 



V3' 

We also have: 



(4.18) 



R.R, 



4d2 



-K„ 



(4.19) 



with 




(4.20) 



the special relation linking the R,'s and vq. 



|Ri + |R2 + IR3 = VQ 



(4.21) 
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Concerning the vector positions of the remaining 9 second neighbors, 3 of them are given by 
— Ri, — -R2/ ~^^3 and the other 6 by the linear combinations R4 = Vq2, R5 = V13, R(, = V03 with 

V02 = Ri + R2, Vi3 = -R2 + ^3, Vo3 = Ri+R2 + R3- (4.22) 

rm and rj^ in the Ai and B4 



From this construction, it follows that generic positions = 
sublattices are given by 

: = ffJjRi + m2R2 + m^R^ , 



= fm + V 



(4.23) 



where m = {mi,m2,mj,) is an integer vector and where the shift vector v = rj^ — is one of 
v,'s in jilTb . 

1) Energy dispersion relation 

First notice that as far as the electronic properties are concerned, the figures (a), (b), (c) of (SJl 
are respectively associated with the sp^, sp^ and sp^ hybridizations of the atom orbitals; i.e: 



figures 


hybridization 


example of molecules 


©a) 


sp^ 


acetylene 




2 

sp 


graphene 




sp^ 


diamond 



(4.24) 



In (6l-a) and ^h), the atoms have delocalized pi- electrons that capture the electronic 
properties of the lattice atoms and have the following dispersion relation. 



^su(N) 



\ 



N-l 



N 



+ 2 ^ cos [flk. (a,-A^ 

i<j=o 



(4.25) 



with N = 2, 3. However, in the case of sp^, the atoms have no delocalized pi-electrons; 



they only have strongly correlated sigma- electrons which make the electronic properties of 
systems based on ilsu(4) different from those based on £51,(3) and ilsu(2)- Nevertheless, as 
far as tight binding model idea is concerned, one may consider other applications; one of 
which concerns the following toy model describing a system based on the lattice ilsu(4) with 
dynamical vacancy sites. 
2) Toy model 

This is a lattice QFT on the £su(4) with dynamical particles and vacancies. The initial state of 
the system correspond to the configuration where the sites of the sublattice A^ are occupied 
by particles and those of the sublattice B4 are imoccupied. 



sublattice 


initial configuration 


quantum state 


A 


particles at r,„ 






vacancy at + v 





Then, the particles (vacancies) start to move towards the neighboring sites with movement 
modeled by hops to first nearest neighbors. Let Ar„, and Br,„+v, be the quantum states 
describing the particle at r,„ and the vacancy at r,,, -|- v,- respectively. Let also and Bj^i^^^ be 
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Fig. 13. the lattice ilsu with sublattices >ls„ (in blue) and S5„(4) (in red). Each atom has 
4 first nearest neighbors, forming a tetrahedron, and 12 second nearest ones. 

the corresponding creation and annihilation operators. The hamiltonian describing the hops 
of the vacancy /particle to the first nearest neighbors is given by 




(4.27) 



By performing the Fourier transform of the wave functions Af^^^, B*^^ , we end with the 
dispersion energy ±t le^l where 

£k = /4 + 2£ cos (k.V,^) , (4.28) 

and V,^ are as in | |4.21I4.22> . The Dirac points are located at ks = ±^ w, with s = 1,2. 
5. Four dimensional graphene 

The so called four dimensional graphene is a QFT model that lives on the 4D hyperdiamond; 
it has links with lattice quantum chromodynamics (QCD) to be discussed in next section. In 
this section, we first study the 4D hyperdiamond; then we use the results of previous section 
to give some physical properties of 4D graphene. 

5.1 Four dimensional hyperdiamond 

Like in the case of 2D honeycomb, the 4D hyperdiamond may be defined by the superposition 
of two sublattices Ai and Bn with the following properties: 

• sites in Aii and S4 are parameterized by the typical M- vectors rn with n = {ni,n2, M3, W4) 
and rij's arbitrary integers. These lattice vectors are expanded as follows 

Ai. In = «i ai + W2 a2 + «3 as + "4 a4 / rj, = rn + 65 , (5.1) 

where ai, a2, as, a4 are primitive vectors generating these sublattices; and 65 is a shift 
vector which we describe below. 



Graphene and Cousin Systems 



23 



• the vector 65 is a global vector taking the same value V n; it is a shift vector giving the 
relative positions of the B4 sites with respect to the Ai ones, 65= — tn, V n. 

The a;'s and 65 vectors can be chosen as 

ai = ei - 65, a2 = 62 - 65, 33 = 63 - 65, 34 = 64 - 65 (5.2) 

with 

= 3 ('+^/5,+^/5,+^/5,+l') , 
4 = \ ( + V^,-V^,-V5,+l) , 

^ * V / (-53-) 

e^' = U-V^,-V^,+V5,+l\ 

4 = \ (-\/5,+\/5,-v^,+lj , 

and 65 = — Ef^i e,-. Notice also that the 5 vectors Cj, 62, 63, 64, 65 define the first nearest 
neighbors to (0, 0, 0, 0) and satisfy the constraint relations, 

e,.e, = 1, e,.e^ = cos !9,y = - \, i ^ ] , (5.4) 

showing that the e,'s are distributed in a symmetric way since all the angles satisfy cos = 
^ ; see also figure fT4t for illustration. 




Fig. 14. On left the 5 first nearest neighbors in the pristine 4D hyperdiamond with the 
properties ||e, || = 1 and ej + £2 + £3 + 64 + ^5 = 0. On right, the 3 first nearest in pristine 
2D graphene with ||e, || = 1 and + 62 + 63 = 0. 

some specific properties 

From the figure fMt representing the first nearest neighbors in the 4D hyperdiamond and 
their analog in 2D graphene, we learn that each An^ type node at rn, with some attached wave 
function A^^, has the following closed neighbors: 

• 5 first nearest neighbors belonging to B4 with wave functions B].^_|_rfe.; they are given by: 



lattice position 


attached wave 


rn + rfei 




In + de2 


Srn+de2 


tn + des 




tn + dCi 




tn + des 


^r„+rfe5 



(5.5) 
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Using this configuration, the typical tight binding hamiltonian describing the couplings 
between the first nearest neighbors reads as 



fn 1 = 1 



(5.6) 



where t is the hop energy and where d is the lattice parameter. 

Notice that in the case where the wave functions at and r,, + de, are rather given by two 
component Weyl spinors 



A" 



Tn+dei 



(5.7) 



together with their adjoints A^^ and , as in the example of 4D lattice QCD to be 

described in next section, the corresponding tight binding model would be. 



5 



4 

Eef 



A" rr'' 



+ }lC 



(5.i 



Tn !=1 L;i = l 

where the e' 's are as in f5.3t ; and where c^, cr^, cr^ are the Pauli matrices and tj^ — ^2x2- 



Notice moreover that the term Yfi=i c,- \A"n^'^a^tn) vanishes identically due to Yji=i = 
0. 

20 second nearest neighbors belonging to the same A/^ with the wave functions 



rn+d(ei-ey) 



they read as 



rn ± d (ei - 62) , fn ± d (ei - 63) , rn ± d (ei - 64) , 
rn ± d (ei - 65) , rn ± d (e2 - 63) , rn ± d (e2 - 64) , 
rn ± d (62 - 65) , rn ± d (63 - 64) , rn ± d (63 - 65) , 
rn ± d (64 - 65) . 



(5.9) 



The 5 vectors ei , 62, 63, 64, 65 are, up to a normalization factor namely , precisely the weight 
vectors Aq, \\, A2, A3, A4 of the 5-dimensional representation of SU (5); and the 20 vectors 
^e, — Cy j are, up to a scale factor their roots /3,y = ^A,- — Ay j . We show as well that the 

particular property e,.ey = — |, which is constant V e, and ey, has a natural interpretation in 
terms of the Cartan matrix of S(J (5) . 
2D/4D Correspondence 

First notice that a generic bond vector e, in the hyperdiamond links two sites in the same unit 
cell of the 4D lattice as shown on the typical coupling term ^r„B+^^^ . This property is quite 
similar to the action of the usual 7'' matrices on 4D (Euclidean) space time spinors which links 
the components of spinors. 

Mim icking the tight binding model of 2D graphene, it has been proposed in l lBedaque et all 
a graphene inspired model for 4D lattice QCD. There, the construction relies on the use 
of the following: 
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the naive correspondence between the bond vectors e, and the 7' matrices 

e, < — > 7/ , = l/-v4 , 

together with 

-65 = 61+62 + 63 + 64 , 

-rs = 7i + 72 + 73 + 74 ■ 

as in the case of 2D graphene, .44-type sites are occupied by left cp'^ and 
2-component Weyl spinors. B4-type sites are occupied by right Xr+de ^^^^ r 
spinors. 



lattice 


2D graphene 


4D hyperdiamond 


^4-sites at r„ 


Ar 




B4-sites at r„ + dej 


R + 




couplings 




< iryLx^de,) 



(5.10) 

(5.11) 

right 
-rfe.Weyl 

(5.12) 



where the indices « = 1, 2 and a = 1,2; and where summation over // is in the Euclidean 
sense. 



For later use, it is interesting to notice the two following: 

(a) in 2D graphene, the wave functions and B^^^^^ describe polarized electrons in first 
nearest sites of the 2D honeycomb. As the spin up and spin down components of the 
electrons contribute equally, the effect of spin couplings in 2D graphene is ignored. 

(b) in the 4D hyperdiamond, we have 4+4 wave functions at each At-iype site or B4-type 
one. These wave functions are given by: 

(i) ip" = ^(fj , j and cp"^ = ^^j: , j having respectively positive and negative 7^ 
chirality, 

<ii> K+de, = {xl+de,' Xr+de) ^^"^ X^+de, = {xl+de/ Xr+de) having respectively negative 
and positive 7^ chirality. 

By mimicking the 2D graphene study, we have the couplings 



u u 

6 ■ a ■ 

I 11 




(5.13) 



building the hamiltonian 

To describe 4D lattice fermions, one considers 4D space time Dirac spinors together with the 
following 7'' matrices realizations, 

-yl = T-*- ® (jl , 7^ = T-*- (g) (J^ , 7^ = T-*- (g) Cr-' , 



7* = T'^ (g) ^2 / 7 
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where the t"s are the Pauli matrices acting on the sublattice structure of the hyperdiamond 
lattice, 

;;)--(roO-^^(;-0' 

The 2x2 matrices a' satisfy as well the Clifford algebra a'af + a) a' = 2S')l2 and act through 
the coupling of left/right Weyl spinors at neighboring sites 

Xl^^,_ - A = {'Pr-%^,f,, - (5.16) 

where at' = (a^ ,(t^, cr" , +il2) and at' = (cr^, a"^, , —ih) ■ For later use, it is interesting to set 
and similar relations for the other cej and a.ej. 

Now extending the tight binding model of 2D graphene to the 4D hyperdiamond; and using 
the weight vectors A, instead of e,, we can build a free fermion action on the 4D lattice by 
attaching a two-component left-handed spinor cp" (r) and right-handed spinor cpl to each 
A4-node r, and a right-handed spinor x" ^ys^ ^^'^ left-handed spinor x" ^-js^, *° every 

B4-node at r + d-^A;. 

The hamiltonian, describing hopping to first nearest-neighbor sites with equal probabilities in 
all five directions A,, reads as follows:. 



r (=0 



Expanding the various spinorial fields in Fourier sums as J * with k 

standing for a generic wave vector in the reciprocal lattice, we can put the field action H4 into 
the form 



where we have set 



D = E Dje'rf^k.A, ^ ^ I ^ ^j'g,rff k.A, j ^ (5 20) 

/=0 fi=l \/=0 



(5.19) 



with 

D, 



* ^A3 + /A4 a1-/A2 
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and pi — k.A/ = X^^, /c^/Aj'. Similarly we have 



1=0 



,,=1 



;=o 



(5.22) 



5.2 Energy dispersion and zero modes 

To get the dispersion energy relations of the 4 waves components (p]^, (p^, x\i their 
corresponding 4 holes, one has to solve the eigenvalues of the Dirac operator l l5.19t . To that 
purpose, we first write the 4-dimensional wave equation as follows. 



-iD 
it) 



(5.23) 



where (p^ = {<pl^,4'k)' /tk = (/Ck'/fk) Weyl spinors and where the 2x2 matrices D, D 
are as in eqs J5.20|5.22t . Then determine the eigenstates and eigenvalues of the 2x2 Dirac 
operator matrix by solving the following characteristic equation. 



det 



f-EO Dn Di2 
-£ D21 D22 
Dn D21 -E 
VD12D22 -EJ 



(5.24) 



from which one can learn the four dispersion energy eigenvalues Ej (k), E2 (k), E3 (k), E4 (k) 

and therefore their zeros. 

1) computing the energy dispersion 

An interesting way to do these calculations is to act on <5.23t once more by the Dirac operator 
to bring it to the following diagonal form 



( 



DD 
V DD 



Xk 



(5.25) 



Then solve separately the eigenvalues problem of the 2-dimensional equations DDcp^^ = E^cp^ 
and DDXk = E^Xk- To do so, it is useful to set 

u (k) = + ,v{k)=d^+ (5.26) 

with = X]/ a'' exp(!rf^fc.A;). Notice that in the continuous limit, we have id^kf, 

u (k) ^ id^ (ki + ;k2), V (k) ^ id^ (k3 + fk4) . (5.27) 
Substituting J5.26> back into J5.20t and <5.22> , we obtain the following expressions. 



DD = 



\u\^ + \v\^ 
2uv 



,2 I . DD = 



|M|2+|zf 

luv 



luv 



(5.28) 
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By solving the characteristic equations of these 2x2 matrix operators, we get the eigenstates 
^k' with their corresponding eigenvalues E^, 



eigenstates 



Y 2|M||t>|'^k ^ Y 2|M||o|'^k 

_ / uv ^,1 I / »P .,2 
Y 2|»||p|^k ^ Y 2|M||p|^k 



eigenvalues 



= |up + \v\ +2\u\ \v 



iif + \vf - 2 ImI \v\ 



E\ = \u\^ + +2 |m| 
£2 = \uf + \v\^ -2\u\ \v\ 



(5.29) 



By taking square roots of Ej., we obtain 2 positive and 2 negative dispersion energies; these 
are 

(5.30) 



E± = +^(|m| ± El = -^(|h| ± \v\Y 

which correspond respectively to particles and the associated holes. 
2) determining the zeros of E± and EJ. 

From the above energy dispersion relations, one sees that the zero modes are of two kinds: 
eI_ = 0, E^ = 0; and E^ = but E^ = E^^^^ ^ 0. Let us consider the case E^ = E^ = 0; 
in this situation the zero modes are given by those wave vectors Kf solving the constraint 
relations u (Kf ) = v (Kf ) = 0. These constraints can be also put in the form 



Kf.Ai 



(5.31) 



for all values o{ pi = 1,2,3,4, or equivalently like d^KpAj = + 2nNi. Notice that 
setting k = Kp + q with small q = ||q|| and expanding D and D, eq <5.23t gets reduced to the 
following familiar wave equation in Dirac theory 



4 



C^^' 

a-l' 



f^]=E 
Xk 



(5.32) 



6. Graphene and lattice QCD 

In this section, we would like to deepen the connection bet ween 2D grap hene and 4D 
lattice QCD. This connection has been first noticed by M.Creutz jCreutzll200i) and has been 
developed by several authors seen its convenience for numerical simulations in QCD . 

6.1 More on link graphene/lattice QCD 

2D graphene has some remarkable properties that can be used to simulate 4D lattice QCD. 
Besides chirality, one of the interesting properties is the existence of two Dirac points that can 
be interpreted as the light quarks up and down. This follows from the study of the zero modes 
of the 2x2 Dirac operator which corresponds also to solve the vanishing of the following 
energy dispersion relation 

; ; 

^ cos aki + i sinflfc/ = , (6.1) 

1=1 1=1 
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which has two zeros as given by d2.22t . 

To make contact with lattice QCD, we start by recalling the usual 4D hamiltonian density of a 
free Dirac fermion Y = [xp^, ip^, XirXi) living in a euclidian space time. 




where 7'' are the usual 4x4 Dirac matrices given by l |5.14t . Then, we discretize this energy 
density H by thinking about the spinorial waves Y (x^, x^) as Yr„ living at the r„-nodes of a 

four dimensional lattice L4 and its space time gradient like | {^T„+aj{ — ^r,,)- The field 
-^ap is the value of the Dirac spinor at the lattice position r„ + aji with the unit vectors fi 
giving the four relative positions of the first nearest neighbors of r„ . Putting this discretization 
back into ^6^, we end with the free fermion model 

H=J-aL(^L [Yr„7''Yr„+n^ " Y,„+„^7'' Yr„] j ■ (6.3) 

The extra two term Yr„rYr„ and (Yr„rYr„)^ with T = | E^, 7'' cancel each other because of 
antisymmetry of the spinors. Clearly, this hamiltonian looks like the tight binding hamiltonian 
describing the electronic properties of the 2D graphene; so one expects several similarities for 
the two systems. 

Mapping the hamiltonian i6.3) to the Fourier space, we get H = X]k (^k^^Yk) with Dirac 
operator 2? = ^ I]^i=i 7'' sin (flfc^) ,where we have set kp = (k.^); giving the wave vector 
component along the ^t-direction. The V- operator is a 4 x 4 matrix that depends on the wave 
vector components {ki,k2,k3,k4) and has 2* zeros located as 

fci = 0, f ; fc2 = 0, f ; fc3 = 0, f ; /C4 = 0, f . (6.4) 

However, to apply these formalism to 4D lattice QCD, the number of the zero modes of the 
Dirac oper ator should b e two in order to interpret them as the light quarks up and down. 
Following jCreutzl.l200^ . this objective can be achieved by modifying l l6.3t so that the Dirac 
operator takes the form 

■ 4 ■ 4 

23 = 1 £ ^Z- Sin {akp) + ^ y'f cos (flfc,,) (6.5) 
" ;i=i " f'=i 

where j'f is some 4x4 matrix that is introduced in next subsection. 
6.2 BoriQi- Creutz fermion s 

Following dCapitani et al.l 12009 allbh and using the 4-component Dirac spinors Yr = 

the Borifi-Creutz (BC) lattice action of free fermions reads in the position space, by dropping 

mass term mg, as follows: 

Hbc-jhL\L "frV'^r+af. " E Y,+,^Y?'Y J - f £ Y.FY, (6.6) 



30 Graphene Simulation 

where, for simplicity, we have dropped out gauge interactions; and where Y'' = 7'' + 17'''; 
which is a kind of complexification of the Dirac matrices. 

Moreover, the matrix Y appearing in the last term is a 4 x 4 matrix linked to 7'', 7''' as follows: 

4 

7'/' = r - 7^', 2r = ^ 7'', 7'' + iy'l' = Yt', yf'Y + = 2S^"'^' , (6.7) 

p=i 



(6.8) 



Mapping <6.6t to the reciprocal space, we have 

k 

where the massless Dirac operator Cgc is given by 

^BC = + ( Y,, - Y;, ) COS (flfc,, ) 

+ ^(Y„+Y,)smK,)-|r. ^'-"^ 

Upon using Y^, + Y^, = 27^, and Yj, — Yj, = 227^,, we can put Vgc in the form 

Vbc = Dk + Dk - -r (6.10) 
a 

with 

Dk = i (^p7>'smakf}j , Dk = I £ 7''' cos flfc,, j , (6.11) 

where fc^, = k.^f. In the next subsection, we will derive the explicit expression of these fc^/s in 
terms of the weight vectors A; of the 5-dimensional representation of the SU (5) symmetry as 
well as useful relations. 

The zero modes of ^^e points in the reciprocal space; they are obtained by solving VgQ = 
0; which leads to the following condition 

7?' (sinflK,, - cosaK^) - T [ 2 - X] cosaK^, ) = . (6.12) 

,,=1 V /'=i / 

This condition is a constraint relation on the wave vector components Kj,; it is solved by the 
two following wave vectors: 

pomtKBC-Ki=K2 = K3 = Ki = , 

that are interpreted in lattice QCD as associated with the light quarks up and down. 
Notice that if giving up the 7^,- terms in eqs <6.6l6^ ; i.e 7J, — > 0, the remaining terms in Dg^ 
namely Dk ~ 7'' sinaKf, have 16 zero modes given by the wave components i<C^ = 0, n. By 
switching on the 7j,-terms, 14 zeros are removed. 



6.3 Hyperdiamond model 

The hamiltonian Hgc is somehow very particular; it let suspecting to hide a more fundamental 
property which can be explicitly exhibited by using hidden symmetries. To that purpose. 
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notice that the price to pay for getting a Dirac operator with two zero modes is the 
involvement of the complexified Dirac matrices V', Y'' as well as the particular matrix Y. 
Despite that it violates explicitly the SO (4) Lorentz symmetry since it can be written as 



,<=i 



1 
1 

VI/ 



(6.14) 



the matrix T plays an important role in studying the zero modes. The expression of the matrix 
r {6.71 should be thought of as associated precisely with the solution of the constraint relation 
2T — 7'' = that is required by a hidden symmetry of the BC model. This invariance is 
precisely the SU (5) symmetry of the 4D hyperdiamond to be identified below. Moreover, the 
BC hamiltonian Hg^ lives on a 4D lattice L^*- generated by ;< = v ; i.e the vectors 



Vi = 







4 




(4\ 

y 

^3 






, V2 = 














^3 










W2J 




\4! 







(6.15) 



These //-vectors look somehow ambiguous to be interpreted both by using the analogy with 
4D graphene prototype; and also from the SU (5) symmetry view. Indeed, to each site r 6 L|'' 
there should be 5 first nearest neighbors that are rotated by S!i (5) symmetry. But from the BC 
hamiltonian we learn that the first nearest neighbors to each site r are: 




(6.16) 



The fifth missing one, namely r — > i+av^ may be interpreted in the BC fermions as associated 
with the extra term involving the matrix F. To take into account the five nearest neighbors, we 
have to use the rigorous correspondence f — v^, and — > V5 which can be also written in a 
combined form as follows F'^ — )• with F*^ = (rf',F^) and = (v^,, V5) . Because of the 
SU (5) symmetry properties, we also have to require the condition vj + V2 + V3 + V4 + V5 = 
characterizing the 5 first nearest neighbors. To determine the explicit expressions of the 
matrices F^ in terms of the usual Dirac ones, we modify the BC model l |6.6t as follows 



/ 5 5 

"BC Tit 2-1 \ 2-1 "r+flVM ~ 2-1 "r+avM^ ~r 

r \M=l M=l 



(6.17) 



exhibiting both SO (4) and SU (5) symmetries and leading to the following free Dirac operator 



^ = E (r?' + r?,) sin (flfc,,) + ^ (Fs + Fg) sin (flfcg) 



4 



(6.18) 



h E i^i' ~ T/O COS (flfc,,) + 2^ (Fg - Fg) cos (flfcg) 
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where fc^vi = k.VM and where n|^^^ ^g'fl'^M j — \^ Em=i = expressing the conservation 
of total momenta at each lattice site. Equating with <6.9I6.10I6.1H , we get the identities 

Y,, + Y,, = r,, + I^„ Y,, - Y,, = Tf, - T,„ (6.19) 

and 

^ (Fs + fg) Sin (fl/cs) + ^ (Fg - fs) COS (aks) = -%T. (6.20) 
Eqs <6.19t are solved by F^, = Y^,; that is F^, = 7?* + ; (F — 7'') while 

F5 = -2iT for sin (afcg) = , 

F5 = -2F for sin(flA:5) = 1 . ^ ' 

where fcg = — {ki +k2 + k3 +ki). In this 5-dimensional approach, the ambiguity in dealing 
with the /(-vectors is overcome; and the underlying SO (4) and SU (5) symmetries of the 
model in reciprocal space are explicitly exhibited. 



7. Conclusion and comments 

Being a simple lattice-carbon based structure with delocalized electrons, graphene has been 
shown to exhibit several exotic physical properties and chemical reactions leading to the 
synthesis of graphene type derivatives such as graphAne and graphOne. In this book chapter, 
we have shown that graphene has also very remarkable hidden symmetries that capture basic 
physical properties; one of these symmetries is the well known SU (2) invariance of the unit 
cells that plays a crucial role in the study of the electronic properties using first principle 
calculations. Another remarkable hidden invariance, which has been developed in this work, 
is the SU(3) symmetry that captures both crystallographic and physical properties of the 
graphene. For instance, first nearest neighbors form 3-dimensional representations of SLf (3); 
and the second nearest neighbor ones transform in its adjoint. Moreover, basic constraint 
relations like vi + V2 + V3 = is precisely a SLf (3) group property; and its solutions are 
exactly given by group theory. Furthermore, the location of the Dirac zero modes of graphene 
is also captured by SU (3) seen that these points are given by ±|jai, ±|ja2/ ^^"^3 where 
the a,'s are the SU (3) roots that generate the reciprocal space. 

On the other hand, from SU (3) group theory's point of view, graphene has cousin systems 
with generic SU (N) symmetries where the integer N takes the values 2,3,4, .... The leading 
graphene cousin systems are linear molecules with hidden SU(2) invariance; this is precisely 
the case of poly-acetylene, cumulene and poly-yne studied in section 4. The graphene 
cousin systems with hidden SU(4) and SU(5) symmetries are given by 3D diamond; and 
4D hyperdiamond which has an application in 4D- lattice QCD. 

Finally, it is worth to mention that the peculiar and unique properties of graphene are expected 
to open new areas of applications due to its important electronic, spintronic, mechanical 
and optical properties. The challenge is find low-cost-processes for producing graphene and 
graphene-based structures and to tune its properties to the targeted applications such as the 
replacement of silicon in the field of new-type of semiconductors and new electronics, new 
data-storage devices, new materials with exceptional mechanical properties and so on. 
Various attempts are also made to incorporate other atoms within the structure of graphene 
or combine the graphene-based structures with other materials in sandwich type structure or 
in chemical way by binding it to various molecules with divers topologies and fimctionalities. 
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